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Abstract 

^ : 

The two-dimensional inhomogeneous zeta-function series (with homogeneous part of the 

CD : 

^ most general Epstein type): 

• • ■ 

> ■ 

•i-H , 

^ (am + bmn + cn + q)~ s , 

. rn,n<=/j 

is analytically continued in the variable s by using zeta-function techniques. A simple for- 
mula is obtained, which extends the Chowla-Selberg formula to inhomogeneous Epstein 
zeta-functions. The new expression is then applied to solve the problem of computing the 
determinant of the basic differential operator that appears in an attempt at quantizing grav- 
ity by using the Wheeler-De Witt equation in 2+1 dimensional spacetime with the torus 
topology. 
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1 Introduction 



In a recent publication |IJ, dealing with the approach to (2+l)-dimensional quantum gravity 
which consists in making direct use of the Wheeler-De Witt equation, Carlip has come 
across a rather involved mathematical problem. Of course, none of the aproaches that have 
been employed for the quantization of gravity is simple, for different reasons. Here we will 
concentrate only in the specific point of the whole problem that has been risen by Carlip, 
and which concerns the calculation of the basic determinant that appears in his method for 
the case of the torus topology. It is the determinant corresponding to a differential operator, 
D , which has the following set of eigenf unctions and eigenvalues (for explicit details, see M 
and references therein): 

4-7T 2 

\mn >= e 2ni{mx+ny \ l mn = \n - mr\ 2 + V , (1) 

where m and n are integers, and r and r 2 are the usual labels corresponding to the standard 
two-dimensional metric for the torus 



with x and y angular coordinates of period 1 and r = T\ + IT2 the modulus (a complex 
parameter [|lj). Vq is the spatial integral of the relevant potential function [|J. 

At that point, the physical difficulty has boiled down to a well formulated mathematical 
problem which, unfortunately, has no straightforward solution from, e.g. the zeta functions 
which commonly appear in physical or mathematical references. Since let us recall that, in 
fact, the best way to obtain the determinant once the spectrum of the operator is known 
is through the calculation of the corresponding zeta function, (d . After simplifying the 
notation a little, one easily recognizes that one has to deal here with a series of the form 



the prime meaning that the term with both m = n = is absent from the sum. Of course 
this distinction needs not to be done when q 7^ (the value of such term being then trivially 
q~ s ), but it is certainly important for considering the particular case q = (see later). One 
is interested in obtaining the function F(l; a, b, c; q) of a, b, c, q, since this expression comes 
inside a functional integral which involves the relevant variables of the problem. As it stands, 
Eq. ([|) has no sense for s = 1, and it is also clear that analytic continuation to such value of 



ds 2 = t 2 1 \dx + rdy\ 2 , 



(2) 



F(s; a, b, c;q)= (am 2 + bran + cn 2 + q) 



(3) 
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s hits a pole and, therefore, must be conveniently defined. This has been done successfully 
in the literature (see 0). 

In what follows we will calculate the sum (^) — and its corresponding analytical continuat- 
ion — and also its derivative with respect to s, by means of zeta-function techniques. The 
final expresion will be remarkably simple, involving just (apart from finite sums) a quickly 
convergent series of exponentially vanishing integrals (of Bessel function type). It is, in fact, 
a generalization of the celebrated (by the mathematicians) Chowla-Selberg formula. We will 
proceed step by step, starting from some particular, more simple cases. 

The case 6 = will be treated in Sect. 2, and the general homogeneous case (q = 0) 
in Sect. 3 — this is the Chowla-Selberg formula itself. In Sect. 4 we will derive the new 
formula, which is capable of dealing with the general situation (q ^ 0). The explicit use 
of the formula in the quantization of (2+l)-dimensional gravity will be treated in Sect. 5, 
where specific results for this physical application will be given. For the sake of comparison, 
we present in Sect. 6 an alternative treatement by means of Eisenstein series. Finally, Sect. 
7 is devoted to conclusions. 



2 Case b=0 

Remember that 

OO 1 OO 

C00 = E- Cff(s,a) = E ? T w a^0,-l,-2,..., Re s > 1, (4) 

n=l n n=0 \ n i" a ) 

define the Riemann and Hurwitz zeta functions, respectively, which can be analytically 
continued in s as meromorphic functions to the whole complex s-plane (with just a simple 
pole at s = 1). In general, in zeta-function regularization we are done when the result can 
be expressed in terms of these simple functions. However, in the more difficult case involved 
here, we will have a different function as the most elementary one (see Eq. ([H]) below). 

Considering again the general series @, the parenthesis in this expression must be visu- 
alized as an inhomogeneous quadratic form: 

Q(x,y) + q, Q(x,y) = ax 2 + bxy + cy 2 , (5) 

restricted to the integers. We shall assume all the time that b > 0, that the discriminant 

A = 4ac - b 2 > 0, (6) 
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and that q is such that Q(m, n)+g^O,Vm,nG Z. In terms of the corresponding physical 
constants, as they appear in [0, these conditions are indeed satisfied (for the physically 
relevant cases). We start by studying some particular situations which, for the benefit of the 
reader, may be interesting to recall. 

The case 6 = corresponds to a situation that we have considered in former papers (even 
in more general terms) and for which we have already derived explicit formulas @] (see also 
H). In particular, for the series 

oo _ 

ai, a 2 ; ci, c 2 ) = J2 [«i(^i + ci) 2 + a 2 (n 2 + c 2 ) 2 + c . (7) 

ni,n.2=0 

we have obtained the following expression: 

a 2 S ^ (-l)"T(s + m) MA™ , 
E 2 {s; ai, a 2 ; c u c 2 ) = -— ^ — — Ctf(-2m, a) 

1 \ S ) m=0 m - VCt2/ 

xEt /a2 (s + m; 1; c 2 ) + ^-1 J^ T ^(s - 1/2; 1; c 2 ) 



O S oo oo 

. ^ —a/2— 1/4 -s/2+1/4 s- 1 / 2 /q > 

+ p7~T a i a 2 n i cos(27rniCiJ 

1 V s / m=l n 2 =0 

27TOi 



a 2 (n 2 + c 2 ) 2 + c 



-s/2+1/4 



/ lb l / o \ 

Xi^s-1/2 I — ^V" 2 ^ 2 + C2 ) +c ) ' ( 8 ) 



where 



rr s ^ (-l) m T(s + m) fai\ m , , n v c 1 / 2 - 5 /¥" T (s - § 



^(s;ai;ci) = — — £ ^ — ^ ^ " — Cff(-2m,ci) + 



r(s) ^ m! V c J 2 y ai r(s) 

+ 7^ ar /2 " 1/4 c- s / 2+1 / 4 £ n s r V2 cos(2vrn lCl )^s-i/2 (^nu^) , (9) 

and being K v the modified Bessel functions of the second kind. In order to specify the 
formula por the present case, we just have to substitute: c\ = c 2 = 1 (the Hurwitz zeta 
functions turn simply into ordinary Riemann ones), ai = a, a 2 = c and c = q, and be a bit 
careful with the summation range; but this can be easily taken care of (see, for instance, 
M where the explicit formulas relating the cases of doubly infinite summation ranges and 
simply infinite ranges are given). 



3 Case q=0 

Here, the Chowla-Selberg formula |J for the (general homogeneous) Epstein zeta function 
|| corresponding to the quadratic form Q is to be used. (This is an expression well known 
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in number theory but not so much in mathematical physics.) The result is 

2 s+3 /V 



F(s;a,b,c;0) = 2((2s) a~ s + 



2 2 \/^a s - 1 



T(s)A s - 1 / 



T( S - 1/2)C(2 S - 1) + 



r( s )A s / 2 -V4^ 



x 



^ n s_1/2 (Ti_2s(^) cos(nwb/a) / dtt s ~ 3/2 exp 



n=0 



7rnv A 
2a~ 



(t + r 1 ) 



(10) 



where 

c£|n 

namely the sum over the s-powers of the divisors of n. 

This formula is very useful and its practical application quite simple. In fact, the two first 
terms are just nice, while the last one (impressive in appearence) is very quickly convergent 
and thus absolutely harmless in practice: only a few first terms of the series need to be 
calculated, even if one needs exceptionally good accuracy. One should also notice that the 
pole of F at s = 1 appears through ((2s — 1) in the second term, while for s = 1/2, the 
apparent singularities of the first and second terms cancel each other and no pole is formed. 

A closer, quantitative idea about the integral can be got from the following closed ex- 
pression for it: 



(12) 



K v being again the modified Bessel function of the second kind. In particular, by calling the 
integral 



I(n,s) 



i 

dtt s ~ 3/2 exp 



■nny A 
2a~ 



(t + t- L ) 



(13) 



one has 



l(n,0) 



Kn,2) 



J(n,3) 



/ 2a 



exp 



7m\/A N 



a + -rcn^/A I 2a 



exp 



= I(n, 1), I(n, 1/2) = 2K (7mVA/a) 

7TOa/A\ 



TrnvA V n\^K 
3a 2 + 3imay/A + n 2 n 2 A / 2a 



TT 2 n 2 A 



n\fA 



exp 



7m-\/A N 



(14) 



As functions of n, all these expressions share the common feature of being exponentially 
decreasing with n. 
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4 The general case a,b,c,q ^ 



This case is more difficult. To handle it, we can choose to go through the whole derivation of 
the Chowla-Selberg formula for the quadratic form Q and see the differences introduced by 
the inhomogeneity (the constant q). Instead, we will here undertake a more down-to-earth 
derivation, which will be similar to the ones that we have successfully employed several times 
in former papers — in particular to obtain Eqs. @ and ([|). Since the technicalities of the 
method have been abundantly discussed before |§ (see also ||), we will here consider the 
main steps of the proof only. They are the following, (i) Rewrite the initial expression (D), 
(|8|), by using the gamma function identity 

£(Q + gr = — -W duu s - l e-^ u . (15) 

m,n *■ \$ ) m,n ® 

(ii) Expand the exponential in terms of power series of m and n and interchange the order 
of the summations, i.e., the sum over such expansion with the sums over m and n or - 
equivalently in this case — use Jacobi's theta function fundamental identity (as can be found, 
for instance, in ||). The equivalence of both methods was explicitly proven in ||. The second 
one starts here from a trivial rewriting of the non-negative quadratic form Q(m,n) as the 
sum of two squares 



Q(m, n) 



' bn\ 2 A , 
m H H -n 



(16) 



2a J Aa? 

and proceeds by considering the summation over m, while treating first n as a parameter, 
(iii) Finally, make the following change of variables (for convenience) 

2irm ,„ . 

u = -=-t (17) 
VA n 



and use the same idea as in Eq. (pTOj) of rewriting the double sum as a sum over the product 
mn and (a finite one) over the divisors of the product: 

E ( Jr) S " /2 = E (n in2 r 1/2 n\-^ = E -*- 1/2 E d l - 2s (18) 

ni,Ti2 ^"2/ n\,n% n d\n 

(this factor appears when the change of variables is performed). On the other hand, the term 
£ n = bn/(2a) in the first square of the decomposition ( |TED (that may be written (m+£ n ) 2 -|-n 2 ) 
leads to a cosine factor in the final expression, e.g. 

cos(27r£ n ) = cos(mib/a). (19) 
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This is also explained in detail in [0 (but notice the small mistake in the first of these 
references, that was later corrected in the subsequent ones). 

By doing all this, the following generalized expression is obtained 

F(s;a,b,c;q) = V] [Q(m } n) + q]~ s = (am 2 + bmn + cn 2 + q)~ s 

= 2Ceh{s, Aaq/ A) a~ s + ^^-1/2 T ( s ~ 1 MCeh(s - 1/2, Aaq/ A) 
2 s+3 /V 



r(5)AV2-l/4^E 
7TO\/A 



n 



;-l/2 



x exp 



2a 



1 + 



cos(mrb/a)Td 1 ~ 2s dtf- 3/2 
« Jo 

Aaq \ 



Ad 2 



d\n 
t + t' 1 



(20) 



where the function (eh{s,p) [one dimensional Epstein- Hurwitz or inhomogeneous Epstein) 
is given by 



(eh{s;p) = J2( n2 +P) 



(21) 



n=l 



v s , >A : r(s-i/2) +1/2 + 2- P 



2T s) 



2ir s r>- s l 2+1 / A 00 

-J2^ 1/2 Ks-i/2(2nn^), 



r s) 



n=l 



and is studied in full detail in [|H| (with numerical tables, plots, and a couple of explicit 
physical applications). 

It is remarkable that the integral inside the series of the new expression can still be 
written in a closed form using (|T2|) — as in the case of Eq. ([H]) . Calling now the integral 

im\f~K 



roc 

J(n,s) = / dtt s ' 3/2 exp 
Jo 

we obtain, in particular, 
J(n, 0) = 



2a 



t + r 1 



(22) 



2a 



V ny/A 
J{n, 1/2) = 2K 



exp 
irn 



Tin ( Aaq\ 

- ( + if) 



i/2- 



'A + 



J(n, 1) 
J(n,2) 



'2aVA 



Aaq 
If J ' 

1/2 



exp 



+ t/A + 



Aaq] / 2aA 3 / 2 



7m 



(i 2 



x exp 



7m 



— A + 



Aaq 



n 
1/2 



a 
A + 



A 



Aaq 



d 2 

4 a g\ -3/2 



1/2' 



7 



J(n, 3) 



3a 2 3a 
+ — 



7T 2 n 2 7TO 



'2aA 5 /2 



x A + 



4ag 
~d? 



-5/2 



exp 



Tin 



— A + 



4ag 



n 

1/2- 



(23) 



which are again exponentially decreasing with n. 

Expression (p4|) itself can be written also in terms of these Bessel functions: 

F(s; a, b, c; q) = 2( EH (s, Aaq/A) a~ s + T^^i, T ( s ~ ^)Ceh(s - 1/2, 4ag/A) (24) 



r(s)A s -V2 



8(2tt) s « a _ m ,/^,i- 2s / A ,%\ 



5/2-1/4 



7m 



4ag 



r( s )A-v^f =0 

Eq. ( [24] ) is the fundamental result of this paper and must be given a name. We propose to 
call it inhomogeneous or generalized Chowla-Selberg formula. To our knowledge, it has never 
appeared before in the mathematical (or physical) literature. 



5 Explicit application of the formula in quantizing grav- 
ity through the Wheeler-De Witt equation 

As discussed in [|IJ, the quantization of gravity in 2+1 dimensions by means of the Wheeler- 
De Witt equation, in a spacetime with the topology RxT 2 (T 2 being the two-dimensional 
torus), of standard metric given by (0), proceeds through the calculation of the zeta function 
corresponding to the basic differential operator Do, which has a spectral decomposition given 
by (0). In terms of the function F(s; a, b, c; q) (§), the zeta function of Dq is 

( Do (s) = F (s; 4tt 2 /t 2 , -8n 2 n/r 2 , 47r 2 (r 2 + r 2 2 )/r 2 ; V Q ) . (25) 

One has, in particular, A = 647T 4 and using Eq. ( f2~4"|) one gets 

2 -2 S +l 7r -2 S 



(d (s) = - (eh UV(47tV 2 



r 2 s 



2 -2 S+ l 7r -2s + l/2 r / s _ 1 / 2) 

+ jzi^ — Ceh (s - 1/2, K,/(47r 2 r 2 ) 

OO />OC 

E^ _1/2 cos(2n 7 rr 1 ) E^~ 2s / ^~ 3/2 



x exp < — niTT2 



i • -4^H 

47r 2 d 2 r 2 



(26) 
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with 



/■ ( t/ ///i 2 \\ 9 2s-l 2s/^o\ " o2s -2 2 S -l/2 r ( S ~ V 2 ) /K)^ 5+1/2 

CBH [S, Vo/i^TT T 2 )) = -2 7T ^ — J +2 



7T 



r s 



^"2 



+ 



23+1/2^23-1/2 /y Q \ -s/2+1/4 oo 



r s 



72 



n=l 



(27) 



The quantity of interest is the determinant of the operator D (see [|IJ). This is most 
conveniently computed by means of its zeta function. In particular: 



det 1/2 £> = exp 



I" 
2 



(28) 



Thus, we must now calculate the derivative of (p4 ) at s = 0. We have, for the general 
function F(s; a, b, c; g) 



F'(0;a,b,c;q) = ]na + 2£ , EH (0,4aq/A) £ EH (-l/2,4aq/A) 



where 



+4 ^2 n 1 cos(niib/ a) ^ d exp 

n=l d|n 



1 



1/2' 



; 



Cbh(0; p) = -vr^ + - lnp + 2p^ ^ n- 1 / 2 K 1/2 (2n7r v ^), 

^ n=l 



(29) 



(30) 



PPC^(-l/2;p) = -| - ^ - ^ £n- 1 K 1 (2n 7 rVp). 



while for (eh(— 1/2; p) the principal part prescription (PP) is to be used (see g, [13], |14|1): 

(31) 

* '« n=l 

Finally, for the determinant of D , we obtain 

det 1/2 /J = exp j-1 lna - ('eh (0| V(4tt 2 t 2 )) + 2ttt 2 C^ (-1/2; y /(47r 2 r 2 

V„ \ 1/2 



— 2 ^ n 1 cos(2n7rri) ^ o? exp 

cJ|tt. 



n=l 



-2nTTT2 1 



47r 2 r 2 d 2 



(32) 



with C'eh (0; K)/(47r 2 r 2 )) and G?h (-1/2; V r /(47r 2 r 2 )) being given by expressions (^) and (0 
above, putting p = Vo/(47r 2 r 2 ). This yields 



JV 



exp < — - — h 



Vb 1 Wo 1 



7T 2 V r 2 2 

9 /T/ \ 1 / 4 00 / 

E-- 1/2 ^i/ 2 

n=l ' 



72 Vb 



_0 

7T \T 2 



— 2 V] ri 1 cos(2n7rri) V] (i exp 

n=l d|n 



T 2 y 77 n=l 



-2n7TT 2 1 + 



TaVbX) 71 ^l I " i 

Vb A 1/2 1 
47r 2 r 2 (i 2 



(33) 
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We observe again that the final formula is really simple since, in practice, it provides a 
very good approximations with just a few terms, which are, on its turn, elementary functions 
of the relevant variables and parameters. This is so, because the infinite series that appear 
converge extremely quickly (terms exponentially decreasing with n). In an asymptotical 
approach to the determinant, only the first line in Eq. (|33|) is relevant (as we will show) and 
the three series can be eliminated altogether. 

From the detailed analysis in Jl|, it follows that the quantity to be calculated now is 
the derivative with respect to Vq of the above determinant, since this quantity vanishes pre- 
cisely at the solutions of the Hamiltonian constraint (always in the language of quantization 
through the corresponding Wheeler-De Witt equation). In other words, the solutions of the 
equation 



_d_ 



det 1/2 L> = 0, 



(34) 



will yield the conditions that the quantized magnitudes and parameters are bound to satisfy 
as a consequence of the Wheeler-De Witt equations. If this does not provide all the solutions 
of such (very involved) differential equations, at least gives us important clues about their 
behaviour (this is one of the two basic problems of the approach in [p], namely that of 
understanding the determinant of Dq). After some calculations one finds that Eq. (^) can 
be written as 



det 1/2 £> 



o x 



1 



57T 



1 

2V 



,3/4 oo 



n=l 



2V27TT 2 
/ 27TT 2 „ | 



(r 2 /V c ^ 4 ~ 



Vq 



/2 



1 

4 
n\ 



1 



— E^i \ n \ - 



n=l 



T 2 



(35) 



1 

27 



cos(2r?,7rr 1 ) E ( d 



n=l 



d\n 



V 



4ir 2 T 2 



-1/2 



exp —2mrT 2 \ 1 



Vn 



An 2 T 2 (p) 







(the primes mean here derivatives of the Bessel functions). In principle — the consistency 
of the approximation is to be checked a posteriori — Eq. ([35]) can be reduced to the very 
simple expression 

1 V f 1\ [Vo 2 

+ 1 W— + -~0. (36) 

27T T 2 V r 2 / V r 2 t 2 
The analysis of this last equation is easy to do. From its discriminant it turns out that real 

solutions can only be obtained when 



r 2 < T. 



(i) 



or 



. (2) 
72 > ?2 ) 



(37) 
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with 

2 2 
r 2 (1) = 1 + - (1 - V^TT) = 0.34104103, r 2 {2) = 1 + - ( 1 + VtF+T) = 2.93219852. (38) 

Moreover, the special situation when r 2 = 1 leads to the non-physical result Vq = —Ait. 



Now, from the constraint y Vq/t 2 > (in order the whole approximation to have sense) we 
see that the second possibility in fl3~7|) just disappears. For the first, when r 2 — > it turns 
out that JVq/t2 ~ 27r/r 2 — > oo or JVq/tz — > 2. When T2 moves within the allowed range 
< r 2 < r^ 1 ^ (the first interval of Q37|)), the corresponding two solutions \JVqJt2 of the 
quadratic equation (1361) sweep the following intervals 



6.07, 



< r 2 < r 2 (1) = 0.34 — > { ff + 2^_2v%+T (3g) 

in the order indicated, respectively. It is easy to check that the first of these two alternatives 
(corresponding to the bigger root of (|36D) provides an absolutely consistent approximate 
solution to the exact equation (^). On the contrary, the second alternative (corresponding 
to the smaller root of (pop) does not actually provide a consistent approximate solution 
(unless r 2 is close to r 2 ). Direct investigation of other possible roots of Eq. (^) (e.g., 
involving some terms of the series) is a quite difficult issue. 



6 An alternative treatement by means of Eisenstein 
series 

An alternative way of treating the general case is the following (see [0]). The inhomogeneity 
(the q term here) is taken care of by the simplest (but hardly economic) method of performing 
a binomial expansion of the sort |14| 



t^T^^E^s + k), (40) 
fc =o k]T ( s ) 



where E(z,s) is an Eisenstein series (see, for instance, Lang |TT| or Kubota |T2(]), which is 
obtained from F(s; a, b, c; 0) by doing the substitution 

u 

2z = a + iu, c = C-, (41) 

2 



so that 



E(z,s)= J2 (u/2) s \m + nz\- 2s , (42) 

m,n=0 
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and has the series expansion 

*(,,.) = 2C(2 8 ) + 2W( M /2)'-' r(1 -^f- 1) 

/2lnlV~ 1/2 

+ 2 E E e ™ a K a _i /2 (mrmu/2). (43) 

m=l n^O V mM / 

It is important to notice, however, that when doing things in this last way the final result 
is expressed in terms of three infinite sums, while in the first general procedure only one 
infinite sum appears (together with a finite sum, for every index n, over the divisors of n), 
and it is very quickly convergent. Notice, moreover, how the d-term in the exponent in Eq. 
(|24|). when expanded in power series, gives rise to the binomial sum corresponding to the 
last treatement. The advantage of the use of the method developed in Sect. 5, steming 
from Eq. (|24l), seems clear (expanding a negative exponential is in general computationally 
disastrous) . 



7 Conclusions 

The main results of this paper have been the derivation of equation (|24]) which generalizes 



the Chowla-Selberg formula ( |iH| ) and its physical application to calculate the determinant 
(|32|) and its derivative, Eqs. (|35|) and (p6[). To have at our disposal an exact expression 
for dealing with inhomogeneous Epstein zeta functions is certainly an interesting thing from 
the physical point of view, since this kind of zeta functions appear frequently in modern 
applications that are no more restricted to zero temperature or Euclidean spacetime. The 
simplicity of the result is remarkable. Namely, the fact that, in practice, we just need a couple 
of terms of the formula ([24]) even if we want to obtain very accurate numerical results. 

Two problems were singularized out in Ref. |l[] as the main difficulties that appear in 
the quantization of (2+l)-dimensional gravity through the Wheeler-De Witt equation: (a) 
to give grounds for the choice of the specific operator ordering of the Hamiltonian constraint 
which leads to the Wheeler-De Witt equation of the quantized system, and (b) to understand 
the functional dependence of the determinant det 1/2 D in terms of the relevant variables 
and to obtain its extrema as a function of the potential Vq. In the present paper we have 
been able to solve problem (b) — at least partially — by means of a numerically consistent 
approximation. Furthermore, we have now also the possibility, through Eqs. (|35| ) and (p6|). 
of getting relevant hints which could lead to the solution of problem (a) — at least in a crude 
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way. In the sense that if, once analyzed in detail, the constraints (|39|) turned out to be 
unphysical, that should lead us to conclude that probably the operator ordering chosen in 
the Hamiltonian constraint (and thus the Wheeler-De Witt equation itself which emerges 
from it) ought to be modified. 

To be noticed is also the more technical point that, in general, when performing the 
analytic continuation through s (necessary, e.g., for the calculation of determinants of dif- 
ferential operators), or when taking the derivative of the series function F with respect to 
s at some particular value of s, it may well happen that we hit a pole since, in general, 
this continuation is a meromorphic function of s. That occurs in our case for s = 1. Such 
a situation can be dealt with in the usual way, by means of the principal-part prescription 
HT3l |14| . In Ref. 0] (see also @), we discuss several explicit examples (appearing in physical 
theories) where the precise manner of doing it is clearly explained, all the way down to the 
numerical results (see also |10| ) . 

Apart from the physical application that we have here considered — directly dealing 
with the generalization of the Chowla-Selberg formula derived above — in general our new 
expression will be certainly useful in physical situations involving massive theories, finite 
temperatures or a chemical potential in a compactified spacetime. This is the physical 
meaning to be attributed to the constant q. On the contrary, for the mathematical uses 
in number theory the consideration of inhomogeneous quadratic forms for defining zeta 
functions does not seem to be specially relevant. This would explain why such formulas are 
not to be found in the mathematical literature, its derivation being, however, anything but 
straightforward. 
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